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The charged vector p mesons in the presence of external magnetic fields at finite temperature 
T and chemical potential p have been investigated in the framework of the Nambu-Jona-Lasinio 
model. We compute the masses of charged p mesons numerically as a function of the magnetic field 
for different values of temperature and chemical potential. The self-energy of the p meson contains 
the quark-loop contribution, i.e. the leading order contribution in 1/N C expansion. The charged p 
meson mass decreases with the magnetic field and drops to zero at a critical magnetic field eB c , which 
means that the charged vector meson condensation, i.e. the electromagnetic superconductor can be 
induced above the critical magnetic field. Surprisingly, it is found that the charged p condensation 
can even survive at high temperature and density. At zero temperature, the critical magnetic field 
just increases slightly with the chemical potential, which indicates that the charged p condensation 
might occur inside compact stars. At zero density, in the temperature range 0.2 — 0.5 GeV, the 
critical magnetic field for charged p condensation is in the range of 0.2 — 0.6 GeV 2 , which indicates 
that the high temperature electromagnetic superconductor could be created at LHC. 

PACS numbers: 13.40.-f, 25.75.-q, ll.10.Wx 


I. INTRODUCTION 

Quantum chromodynamics(QCD) is widely accepted as the fundamental theory for strong interaction. Exploring 
QCD phase structure and properties of QCD matter at extreme conditions is one of the most important topics of 
high energy nuclear physics. Recently, it becomes more and more attractive to study properties of QCD matter 
under strong magnetic fields. Strong magnetic fields could exist in various physical systems, e.g. the strength of 
magnetic fields can reach up to 10 14 G in magnetars, and the magnitude of 10 18 ~ 23 G of magnetic fields may appear 
in early universe. Especially, in the laboratory, strong magnetic fields with the strength of 10 18 ~ 2 ° G (amounts to 
eB ~ 0.1 — 1 GeV 2 ) can be generated in the early stage of non-central heavy ion collisions 0,1 at the Relativistic 
Heavy Ion Collider (RHIC) or the LHC. Therefore, heavy ion collisions offer a unique platform to probe properties of 
QCD vacuum and hot/dense quark matter under strong magnetic fields. 

Magnetic fields can modify the spectrum of charged particles, the energy levels of free charged particles in an 
external magnetic field parallel to 2 axis are given as E 2 s ^(p z ) = P 2 Z + (2ra — 2 sgn(q)s z + 1)|<Z-B| + m 2 , where the 
nonnegative integer n is the Landau level, s z is the spin projection, and p z is the particle’s momentum in the direction 
of the magnetic field. Hence, m 2 ±(B) = m 2 ±(B = 0) + eB and m 2 ±(B) = m 2 ±{B = 0) — eB describe the minimal 
effective squared masses of charged n and charged p, respectively. As the magnetic field increases, the charged 
7T mesons become heavier while the polarized charged p mesons become lighter. Chernodub firstly proposed that 
charged p meson can form a superconductor state when magnetic field is stronger than the critical magnetic field 
eB c « m 2 p± (B = 0) » 0.6 GeV 2 0 , 0 . This characterizes a quantum phase transition caused only by the magnetic 
field, which is a surprising result. 

It is still under investigation whether there exists a vacuum superconductor. In Ref.0, it was argued that the 
vacuum superconductor phase cannot appear due to the Vafa-Witten theorem. However, the authors in Ref. 0] 
argued that the Vafa-Witten theorem would not forbid the charged p meson condensation in an external magnetic 
field. More studies on the charged p condensation have been done in different frameworks such as the Nambu-Jona- 
Lasinio (NJL) model 0, lattice QCD 0,0], the gauge/gravity correspondence 0,0 , the relativistic Hamiltonian 
technique [12j |. and the Dyson-Schwinger equations [131] . On one hand, it was observed that masses of charged p mesons 
decrease firstly and then increase with the magnetic field, and cannot drop to zero, as shown by lattice calculations 
0, Dyson-Schwinger equations [13] and by solving the meson spectra in a relativistic quark-antiquark system with 


"Electronic address: haoliu@mail.ihep.ac.cn 
^Electronic address: yulang@mail.ihep.ac.cn 
t Electronic address: huangm@mail.ihep.ac.cn 



2 


the relativistic Hamiltonian technique in Ref. [12] ■ On the other hand, some other results showed that the charged p 
masses indeed decrease with magnetic field and becomes massless at a critical magnetic field as shown in the NJL 
model and in the holographic QCD model B3- 

Furthermore, for charged p condensation, the obtained results for the critical magnetic field are quite different in 
different calculations. For example, the polarized p ± condensstion is expected to occur at the critical magnetic field 
eB c ss 1 GeV 2 in lattice QCD [§], and eB c ss 1.08 m 2 (B = 0) is given in Ref. 0 in a holographic approach. Even 
in the framework of the NJL model, the results for the critical magnetic field differ a lot. The critical magnetic field 
eB c > 1 GeV 2 was obtained in Ref.Q by neglecting the quark-loop contribution, and eB c = 0.98 M 2 was given in 
Ref. Q, where M q is the quark mass. In Ref.Q, we have obtained the critical magnetic field eB c ss 0.2 GeV 2 in the 
NJL model by calculating the p meson polarization function to the leading order of 1/N C expansion. 

By taking the quark propagator in the Ritus form as well as the Landau level representation, we have obtained the 
same results for the critical magnetic field eB c ss 0.2 GeV 2 , which is only 1/3 of the results from the point-particle 
calculation in Ref.[3|. Our results suggests that the charged p condensation can be realized in nature much easier. In 
this work, we are going to study the charged p mesons in a magnetic field at finite temperature and chemical potential 
in the framework of the NJL model, and investigate the temperature and density region for the survival of charged 
p condensation. This paper is organized as follows: In Sec. II, we give a general description of the two-flavor NJL 
model including the effective four-quark interaction in the vector channel, and derive the vector meson mass under 
magnetic fields at finite temperature and chemical potential. We give numerical results and analysis in Sec. Ill and 
finally in Sec. IV the discussion and conclusion is given. 


II. MODEL AND FORMALISM 


A. The two-flavor magnetized NJL model 


In this paper, we investigate the properties of the charged p meson under the magnetic fields with nonzero temper¬ 
ature T and chemical potential p in the framework of the SU(2) NJL model. The NJL modeljSeen as a low-energy 
approximation of QCD, is regarded as an effective theory to study the properties of mesons |14B19| . In our model, 
the Lagrangian density is given by 

C = ip(i J/) — to + p^ip + Gs [(ipip) 2 + {4>i^ 5 Tip) 2 ] 

-g v + tfyvW) 2 ] - \f^ v - (i) 


In the above equation, ip corresponds to the quark field of two light flavors u and d , m = diag(m U: nid) is the current 
quark mass matrix of u and d quarks, r a = (/, r) with r = (t,t ,r 3 ) corresponding to the isospin Pauli matrices, 
and Gs and Gy are the coupling constants with respect to the scalar (pseudoscalar) and the vector (axial-vector) 
channels, respectively. The covariant derivative, D M — iqfA couples quarks to an external magnetic field 

B = (0,0, B) along the positive 2 direction via a background field, for example, A/ Xt = (0,0, Bx , 0). qj = (—1/3, 2/3) 
is defined as the electric charge of the quark field. The field strength tensor is defined as usual by = d^A^, 
with Af xt fixed as above. 

r*' 

The semibosonized Lagrangian which is equivalent to the above Lagrangian is given by 


C sb = ip{x) - to + p'y 0 ) ip(x) -ip(a + ij 5 T ■ n) ip 

( ct 2 +7 ?2) | (V“V aiJ - + A’p l A ati ) B 2 
1C,s + 4GV 2~’ 


( 2 ) 


where the Euler-Lagrange equation of motion for the auxiliary fields lead to the constraints 


a(x) = —2Gs < ip(x)ip(x) >, (3) 

tt(x) = -2 G s < ip^i^Tipix) >, (4) 

V£(x) = 2G V < ip(x)'y f j j T a ip(x) >, (5) 

A“{x) = 2G V < ip{x) r jf 1 r )^T a ip{x) > . (6) 


The quarks obtain the dynamical masses by the quark-antiquark condensation, so the constituent quark mass M of 
u and d is defined by 


M = too — 2 Gs < ipip > 


( 7 ) 




3 


where we assume m u = m d = m q. To determine the dynamical quark mass, in other words, the a condensation, we 
need to minimize the effective potential. The one-loop level effective potential of the model is given by 


n = 


4 G s 




Wf eB I 


+ OO 




\ ' f + °° dp 3 

2 ^ a P J 4^2 

p=0 • y -°° 


{/3£ g + In (l + e - /3(E ‘> +/i) ) + In (l + e -(K E «-^ J , 


( 8 ) 

where /? = ^ and a p = 2 — 6 P o is the spin degeneracy factor. 

In the NJL model, mesons are constructed by the infinite sum of quark-loop chains by using random phase approx¬ 
imation. We calculate the p meson polarization function to the leading order of 1 /N c expansion. The propagator 
of the p meson D^(q 2 ) can be obtained from the one quark loop polarization Hp,v, a b{q 2 ) via the Schwinger-Dyson 
equation and takes the form of 

[~iD%] = [-2iGV<WH + [-2 iG v S ac g^] [-iU x ^ cd ] [~iD%] , (9) 

where a, b, c, d are isospin indices and p, v, A, a are Lorentz indices. The one-loop polarization n /iI/ ’ ab (g 2 ) is given 
by 


rp"' ab (g 2 ) = -if d 4 x e ig ' x Tr[7 /i T a 5g(a;, 0)"f v T b S Q (0, x)\. 
Here, Sq{x,d) is the Ritus fermion propagator 


S Q (x,y)=iYpp e-^^Ppix^D^iP) P p ( yi ), 


( 10 ) 


( 11 ) 




arising from the solution of Dirac equation under a uniform magnetic field using Ritus eigenfunction method. In 
Eq.(jn)(, p = (po,0,p 2 ,p 3 ) : Pp = dpo (^)t P3 , and P p (a:i) is given by 


P p(x i) = ^[fp s { x i) + n p/ P s ( a; i)] + l -^r[fp 8 (*i) - n p/ P s (zi)bV 


( 12 ) 


where, n p = 1 — 6 P o considering the spin degeneracy in the LLL, sq = sgn(QeB). The functions f p s {x i) are defined 

by 


fp s (x l) = </>p (xi - s Q p 2 i 2 B ) , p = 0,1,2, 
fp s {x l) = </>p-i (xi - sqp 2 £ 2 b ) , p = 1,2,3, 

where 4> p {x) is a function of Hermite polynomials H p (x) in the form 

<t>p{x) = a p exp (—^ T ) H p . 


2il 


£b 


(13) 


(14) 


Here, a p = (2 p p!-y/7r£g) 1//2 is the normalization factor and £b = \QeB\ 1 / 2 is the magnetic length. In Eg. (fill) . 
Dq{P) = 7 • Pq — M with the Ritus four momentum p = (po, 0, — SQ\j2\QeB\p,p 3 ). Note that Q is a 2 x 2 matrix in 
the flavor space, so the functions f p s (x i) are matrices in the flavor space. 

Using the propagator in Eq. ED, the one quark loop polarization function H /J ,^ a b(q 2 ) is given by 


where 


IP^V) = / VpVk d i xe-^ p -' k - q) - x ^ a \p : k,x 1 ), 

p,k—0 ^ ^ 


A^ ab (p,k, Xl ) = Tr[^T a P p (x 1 )DQ 1 (p)P p (0)^T b K k (0)DQ 1 (k)K k (x 1 )]. 


(15) 


(16) 


We can get the boundary conditions qi = pi — ki (i = 0,2,3) by integrating over Xq,x 2 ,x 3 and then over po,p 2 ,p 3 
components in Eq. dlSl) . So we get 


W v ’ ab (q 2 ) = i 

p,k =0 1 


dk 3 dk 3 

(27r) 3 


dk 2 dxie iqiXl A^’ ab (p, k,x i) 


(17) 


In above equation Eq. ED, the b.c is short for boundary conditions. 
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B. Charged p meson in a magnetic field at nonzero T and p 


In the rest frame of p , the Lorentz and flavor structure of the one quark loop polarization allows for the following 
decomposition 

n Z(q 2 ) = [n;(<P)pr + ni(g 2 )pr + ni(g 2 )£' 1, ' + n2( g 2 )«'*u , ']<y a6 , (is) 

where = (1,0,0,0) is the 4-velocity in the rest frame. Here, we define the spin projection operators 

Pi v = = lforp), (19) 

K = -e?er,(sz = liovp), ( 20 ) 

If v = -b»b v , (s z = 0 for p ), (21) 

where = (0,0, 0,1) corresponds to the external magnetic field direction and e^, e 2 are the right and left-handed 
polarization vectors 


e i = ^(0,l,i,0), 

«2 = ^(0,1,-i,0). 

As a consequence, we can rewrite the propagator of p meson in Eq.© as 

Kbtf) = [Di(q 2 )Pr + D 2 (q 2 )Pr + D 3 (q 2 )L^ + D^)uV]S ab , 

where 


DM 2 ) = 


2 G v 


1 + 2G? v n?' 

At last, we use the following gap equations to determine the masses of p meson with different spin projections 

1 + 2G v M i = 0. 

In Ea. lflGl) . the isospin Pauli matrices are r° = M and r b = for charged p ± meson, and M = -^(r 1 ± it 2 ). In 
the rest frame of p meson, i.e., q^ = (M p ±, 0 ), the one quark loop polarization function takes the form of 


( 22 ) 

(23) 

(24) 

(25) 

(26) 


n;±(g 2 ) = * £ 

p,k —0 *" 


dkodks 


dk 2 dx 1 A^ pk (p,k,xi) 


I b.c' 


(27) 


After calculating the U^±(q 2 ) in the rest frame of p meson (the details are given in [j|), we can get the matrix as 
following: 


= 

p ± 


(° 

0 

0 

0 ^ 


(° 

0 

0 

0 

0 

n 11 

n 12 

0 


0 

a 

ib 

0 

0 

n 21 

n 22 

0 

- 

0 

—ib 

a 

0 

Vo 

0 

0 

H-y 


Vo 

0 

0 

c 


(28) 


Here, we define n 11 = U 22 = a, A 12 = — n 21 = ib and n 33 = c. Note that the calculation in (6j was finished in the 
vacuum. However, in this paper, we are going to investigate the charged p condensation at finite temperature and 
density. Therefore, fco should be replaced by the Matsubara frequency summation, which is given in the Appendix [A] 
The matrix elements n 11 = n 22 , n 12 = — n 21 and n 33 for p + at finite temperature are given as following: 

n 11 = n 22 = 24 ^ f f dk 2 dx\{AA + a + + BA~a~}, 

p,k=0 J ' ^ 

n 12 =-n 21 = 24 ^ f t^-T 2 [ dk 2 dxi{8(-is u A~a + ) + A(—is u A + a~)}, 

P ,k= o"' 

n 33 = 24 ^ / wh [ dk 2 d Xl {VB+p+ +CB~p~}. 


(29) 
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The definitions of A, B , C and T> are in the Appendix lAl and A ± , a , B ± and /3 =l= are represented as follows: 

= \ [/ P + '«(*i)/^ i (*i)±n P n fc /---( a:i )/ fc ---( a : 1 )], 

= \ [/ P + '“(o)/ fc +a -(o)±n f) n fc /-'»(o)/ fc - 8 -(o)], 

B ± = \[U p f^(x 1 )f- s -(x 1 )±H k f+^(x 1 )f^(x 1 )}, 

/3 ± = ^n P /+ s - i (o)/- s “(o) ± n fc /+ s “(o)/ fc - Srf (o)]. (so) 

where {s u ,s d } = (sgn(g„eB),sgn(g d eB)} and {g„,< 3 w} = {2/3,-1/3}. 

Combining the expression in Eo. (11811 and the matrix in Ea. (12811 . we can easily find the relation for charged p ± 
meson 

II? = — ( a + b), 

II2 = b - a, 

n 1 = -C. (31) 

As it lias been discussed in [22j, II 2 should be zero in the rest frame of p mesons, guaranteed by the Ward identity. 


III. NUMERICAL RESULTS AND DISCUSSIONS 


A. Parameters 


For numerical calculations, we use the soft 


cut-off functions in Refs. 


/a 


A 10 


A 10 + k 2 * 5 : 


fleB = 


A 10 


A 10 + (fcf + 2\QeB\k) 5 ’ 


(32) 

(33) 


for zero and nonzero magnetic fields, respectively. At finite magnetic field, we sum up to 20 Landau levels, and 
the results are saturated. Following Refs. @!H, the parameters of our model, namely the coupling constants Gs 
and Gy , the current quark mass mo, and the three-momentum cut-off A, are determined by reproducing the pion 
decay constant / ff , the quark mass M, the mass of n, and the mass of p in the vacuum. We obtain A = 582 MeV, 
G 5 A 2 = 2.388, GyA 2 = 1.73, and mo = 5 MeV by choosing f n = 95 MeV, m, = 140 MeV, M p = 768 MeV, the 
vacuum quark mass M = 458 MeV. 


B. Numerical results for charged p ± and discussions 

In our numerical calculations, we firstly solve the gap equation Eq. 0 to obtain the quark constitute mass M, and 
then solve the gap equations of Ea. CC)l) for charged vector mesons. 

Fig.[l]shows the quark mass dependence of magnetic field eB at several fixed temperatures T and chemical potentials 
/i, and Fig. [2] describes the quark mass dependence of the temperature T at p = 0 for several different magnetic fields 
eB. It is noticed that in the current model, there is no mechanism for inverse magnetic catalysis around the critical 
temperature. Therefore, Fig. [T] and Fig. [ 2 ] only show magnetic catalysis behavior of the quark mass. 

In Ref.[&|], we have observed that there is possible condensation of charged p + (p ~) with s z = +l(s, = —1) in the 
vacuum when the magnetic field is stronger than the critical magnetic field eB c ~ 0.2 GeV 2 . The main goal of this 
work is to investigate the melting of the charged p condensation at finite temperature and chemical potential. 

Fig. [3] shows the mass of charged pA as the function of the magnetic field eB at zero temperature for different chemical 
potentials. It is observed from Fig. [3] that at zero temperature, when p is smaller than the vacuum constituent quark 
mass = 458MeV, the behavior of magnetic field dependence of the charged p ± mass almost does not change 
for different chemical potentials, and the critical magnetic field remains as eB c = 0.2GeV 2 . This indicates that the 
charged p condensation can survive at high baryon density and might occur inside compact stars. 
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FIG. 1: Quark constitute mass M as a function of magnetic field eB with different temperatures T and chemical potentials p. 



FIG. 2: Quark constitute mass M as a function of the temperature T with different magnetic field eB and zero chemical 
potential p = 0. 



eB[GeV 2 ] 


FIG. 3: The mass of charged p + (p ) with s z = +l(s z = —1) in magnetic field with fixed T = 0 MeV with different chemical 
potentials p. 


Fig. [4] shows the mass of charged as the function of the magnetic field eB at zero chemical potential for different 
temperatures. It is seen that for different temperatures, the charged mass decreases with the magnetic field and 
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drops to zero at the critical magnetic field. With the increase of temperature, the critical magnetic field also increases. 
This indicates that the charged p condensation can survive even at high temperatures. 



FIG. 4: The mass of charged p + (p ) with s z = +l(s z = —1) in magnetic field with fixed p = 0 MeV and different temperatures 
T. 


At last, in Fig. [5l we show the critical magnetic field eB c as a function of the temperature for p = 0,0.3,0.46GeV, 
respectively. It can be read that when the temperature increases, the critical magnetic field increases. Furthermore, 
it is seen that the critical magnetic field increases with the chemical potential at fixed temperature. It is worthy 
of mentioning that at zero chemical potential, when the temperature is below T = 250MeV, which is almost the 
critical temperature for the chiral phase transition, we can read that the critical magnetic field does not change so 
much comparing with its value at zero temperature. However, when the temperature is higher than T = 250MeV, 
it is found that the critical magnetic field increases linearly with the temperature, and in the temperature range 
200 — 500 MeV, the critical magnetic field for charged p condensation is in the range of 0.2 — 0.6 GeV 2 , which is just 
located inside the range of the magnetic field generated in the non-central heavy ion collisions of LHC. It means that 
high temperature superconductor could be produced at LHC. 



FIG. 5: The critical magnetic field eB c as a function of temperature T with different chemical potential p . 


IV. CONCLUSIONS 

We have investigated the charged p mesons in an external magnetic field at finite temperature and density by using 
the NJL model. The mesons are constructed by summing up infinite quark-loop chains by using the random phase 
approximation. In this paper, we calculate the p meson polarization tensor to the leading order of 1/N C expansion, 
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i.e. the one quark loop. In this process, the constituent quark mass is solved self-consistently with magnetic field at 
finite temperature and density. It is noticed that in our current framework, there is no inverse magnetic catalysis 
mechanism, so the quark mass increases with magnetic field and only shows the magnetic catalysis effect. 

The mass of charged p meson depending on the magnetic field is calculated at finite temperature and density. It 
is found that at fixed temperature and density, the polarized charged p meson becomes lighter with the increase of 
magnetic fields and drops to massless at the critical magnetic field eB c , which indicates that there will appear charged 
p meson condensation when the strength of the magnetic field is greater than the eB c . Moreover, it is observed that 
the critical magnetic field eB c increases with both the temperature and density. Particularly, at zero density, our 
results show that in the temperature region 200 MeV < T < 500 MeV, the critical magnetic field eB c is in the range 
of 0.2 — 0.6 GeV 2 , which indicates that the high temperature superconductor could be created in the early stage of 
LHC. 

However, we have to mention that in our current framework, there is no inverse magnetic catalysis for the quark 
mass. In the next step, we will investigate how the inverse magnetic catalysis [2f| will affect our results on the charge 
p condensation, especially at high temperatures. 
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Appendix A: Integrals of ko 

We introduce the notations of the integrals of ko for p + meson as follows: 

dko P(u) ■ k(d) - M 2 1 


A = i 
A 


2 t t (Po-ul, P )(ko-Ud,k) 2 

(aT 2 + ~ pf 2 ,u) - k 2 2 d )) + P(2,u)k(2,d) 


= + Hi) 


ih, 


(Al) 


K _. /' dk 0 P(u) ■ k( d ) — M 2 + 2p (2jU )fc ( 2,d) l ( . T , , r ^ 

J to (p§-< P )(fco 2 -<4,) ~2 [lh+lh) 


1 


m 2 + - - 2 


1 + P(2,u) k(2,d)) P(2,u)k(2,d) 


ih 


(A2) 


„ . / dkoPoko + P{2,u)k( 2 ,d) + H - M- 1 , 

C = l / - y -j-n—1771-- = o(“l + lI l) 

J 27r (pi - )(fcg - Uj 2 d k ) 2 


+ 


2ko 


u 


- pf2,u) — k(2,d) ) + P(2,u)k(2,d) 


ih, 


V = i 


dk 0 Po k o ~ P(2,u)k( 2 ,d) + fcf ~ M2 
2 tt {pl-ulA{kl-u 2 ) 


= 


■Hi) 


2 k^, — 


s( 




- k 2 
k (2,d) 


■ P(2,u)k(2,d) 


ih. 


(A3) 


(A4) 
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Here, w 2 p = 2\q u eB\p + pj + M 2 , p {2 , u ) = -s u s/2\q u eB\p and p (u ) = (p 0 , 0, -s u y/2\q u eB\p,p 3 ). We use the similar 
definition for k^)i ^( 2 ,d)> ^d.k- Moreover, I\, I I and I 2 are given by: 


h 

h 

h 


dko 1 

2tt kl - ul >p ’ 

dko 1 

2tt fcg - ’ 

dko 

27r ((fc 0 + q 0 ) 2 


1 


— w 


2 ,p )(^0 - w d,fc)' 


(A5) 

(A 6 ) 

(A7) 


As in Ref.jHj], we replace the integration over ko by Matsubara summation according to the prescription 
/ '') -► * T Em=-oo(- • •) and obtain 


Hi = -( 

*/; = -( 


nf{u u , p - p) + nf(uj U} p + p)~l 

rif(uid,k ~ k) + nf(ujd,k + p) -1 
^^d,k 


), 


), 


(A 8 ) 

(A9) 


.j _ _ n>f(u) d ,k - n) _1_ 

2u d ,k {ud,k + M p +) 2 - W 2 jp 

nf(-u d ,k - d) _1_ 

2uJd,k ( T ) 2 

| Tlf(w U ,p - k) _1_ 

( Mp+ T Ldu,p) 2 ^d^k 

n f(-u UtP ~ Ap_1_ 

( A7p+ Wu,p) 2 ^d,k 

(Aio) 

with 

? V 0 ) = 7 , X x • ( Al1 ) 

1 + er 

Here, the notations are only for p + meson, and the similar ones can be used for p - meson. 
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